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“KPipe” is a proposed experiment which will study muon neutrino disappearance for a sensitive 
test of the Am 2 ~ 1 eV 2 anomalies, possibly indicative of one or more sterile neutrinos. The ex¬ 
periment is to be located at the J-PARC Materials and Life Science Facility’s spallation neutron 
source, which represents the world’s most intense source of charged kaon decay-at-rest monoener¬ 
getic (236 MeV) muon neutrinos. The detector vessel, designed to measure the charged current 
interactions of these neutrinos, will be 3 m in diameter and 120 m long, extending radially at a 
distance of 32 m to 152 m from the source. This design allows a sensitive search for v ' M disap¬ 
pearance associated with currently favored light sterile neutrino models and features the ability to 
reconstruct the neutrino oscillation wave within a single, extended detector. The required detector 
design, technology, and costs are modest. The KPipe measurements will be robust since they depend 
on a known energy neutrino source with low expected backgrounds. Further, since the measure¬ 
ments rely only on the measured rate of detected events as a function of distance, with no required 
knowledge of the initial flux and neutrino interaction cross section, the results will be largely free 
of systematic errors. The experimental sensitivity to oscillations, based on a shape-only analysis of 
the L/E distribution, will extend an order of magnitude beyond present experimental limits in the 
relevant high-Am 2 parameter space. 

PACS numbers: 14.60.Pq, 14.60.St 


I. INTRODUCTION 

A number of experimental anomalies consistent with 
neutrino oscillations at a characteristic mass splitting 
around 1 eV 2 hint at the possibility of an additional neu¬ 
trino. These anomalies fall into two categories: muon-to- 
electron flavor appearance, as observed by the LSND pQ 
and MiniBooNE m ei experiments, and electron fla¬ 
vor disappearance, as observed by reactor II M and 
source EHS] experiments. A favored beyond-Standard- 
Model explanation for these anomalies invokes an addi¬ 
tional number of N sterile neutrinos participating in os¬ 
cillations beyond the three active flavors fTTflfnn . These 
“3+iV models” are able to simultaneously describe the 
existing anomalous observations and those measurements 
which do not claim a signal in the relevant parameter 
space [141122] . although there is tension between both 
neutrino and antineutrino measurements and appearance 
and disappearance measurements. 

Muon neutrinos, for example, must disappear if the ob¬ 
served anomalies are due to oscillations involving a light 
sterile neutrino. The lack of observed ry, disappearance 
is a major source of tension in the global fits. In order 
to understand the importance of disappearance mea¬ 
surements, consider a 3+1 sterile neutrino model, with 
the probability for ry, —> v e appearance given by: 

P{v» -7 v e ) ~ 4,\U^ 4 \ 2 \U e 4 2 sin 2 (1.27Am^L/E) . (1) 

The probability for v e and v t , disappearance are, respec¬ 
tively: 

P(v e -4 v e ) ~ 1 — 4(1— |U e 4| 2 )|C/ e4 | 2 sin 2 (1.27Am^L/£;) 

( 2 ) 


and 

P{y„ -4 v„) ~ 1—4(1—|[/ At4 | 2 )|C7 /i4 | 2 sin 2 (1.27Am 41 L/.E). 

(3) 

In these equations, the elements of the mixing matrix, U , 
set the amplitude of oscillation, while Am 41 establishes 
the oscillation wavelength. Within this 3+1 model, a 
global fit to the world’s data, including all anomalies and 
null results, will simultaneously constrain f7 e4 , (7^4, and 
Am 41 . The range of values that U ^ 4 can take on, and 
therefore the oscillation parameters that govern dis¬ 
appearance, can thus be restricted. The present global 
fit for disappearance places the allowed region just 
outside of current bounds. This motivates the construc¬ 
tion of a fast, low cost, and decisive disappearance 
experiment that can confirm or disallow various models 
for sterile neutrinos, as well as inform a range of future 
proposed experiments [231132] , 

In what follows we describe such an experiment, called 
“KPipe”, that can perform a search for disappearance 
that extends well beyond current limits while still being 
low cost. KPipe will employ a long, liquid scintillator- 
based detector that is oriented radially with respect to 
an intense source of isotropic monoenergetic 236 MeV 
v^s coming from the decay-at-rest of positively charged 
kaons (K+ -+ BR=63.55±0.11% [EH))- As the 

only relevant monoenergetic neutrino that can interact 
via the charged current interaction, a kaon decay-at-rest 
(KDAR) source represents a unique and important 
tool for precision oscillation, cross section, and nuclear 
physics measurements (HU [35]. Since the energy of these 
neutrinos is known, indications of disappearance may 
be seen along the length of the KPipe detector as os¬ 
cillating deviations from the expected 1 /R 2 dependence 
in the rate of ly, charged-current (CC) interactions. A 
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FIG. 1: An aerial view from Google Maps (2015) of the Ma¬ 
terials and Life Science Facility layout with a superimposed 
schematic drawing [22] of the first floor, including the target 
station. The proposed KPipe location (shown with a dotted 
contour) is 32 m from the target station and 102° with respect 
to the incident proton beam direction. The detector extends 
radially outward from the target station. 


measurement of such a deviation over a large range of 
L/E would not only be a clear indication for the exis¬ 
tence of at least one light sterile neutrino, but also begin 
to disambiguate among different sterile neutrino models. 


II. THE KDAR SOURCE AND KPIPE 
DETECTOR DESIGN 

The Materials and Life Science Experimental Facility 
(MLF) at the Japan Proton Accelerator Research Com¬ 
plex (J-PARC) in Tokai, Japan houses a spallation neu¬ 
tron source used for basic research on materials and life 
science, as well as research and development in indus¬ 
trial engineering. It is also an intense, yet completely 
unutilized, source of neutrinos that emits the world’s 
most intense flux of KDAR monoenergetic (236 MeV) 
z/^s. Neutron beams along with neutrinos are generated 
when a mercury target is hit by a pulsed, high inten¬ 
sity proton beam from the J-PARC rapid-cycling syn¬ 
chrotron (RCS) [58] , The RCS delivers a 3 GeV, 25 Hz 
pulsed proton beam, which arrives in two 80 ns buckets 
spaced 540 ns apart. The facility provides users 500 kW 
of protons-on-target (POT) but has demonstrated its 
eventual steady-state goal of 1 MW, albeit for short 
times (3BJ. Along with neutrons, the proton-on-target 
interaction provides an intense source of light mesons, 
including kaons, which usually come to rest in the high- 
A target and surrounding shielding. 

KPipe will search for muon-flavor disappearance with 



FIG. 2: The KPipe detector design, featuring a 3 m inner 
diameter high density polyethylene (HDPE) vessel filled with 
liquid scintillator. Silicon photomultipliers (SiPMs) are seen 
mounted on the interior panels in hoops spaced by 10 cm in 
the longitudinal direction. The cosmic ray veto is a 10 cm 
space between the panels and the outer HDPE wall. 


CC interactions of 236 MeV v^s on carbon nuclei 
(zz M 12 C —> /i _ A') in liquid scintillator. This interaction 
produces a visible muon and X, where A' is some combi¬ 
nation of an excited nucleus, de-excitation photons, and 
one or more ejected nucleons after final state interactions. 
The goal of the KPipe detector design is to efficiently 
identify these 236 MeV CC events, broadly character¬ 
ized by two separated flashes of light in time coming from 
the prompt X followed by the muon’s decay electron. 

The KPipe design calls for a relatively low cost, 3 m in¬ 
ner diameter (ID) steel-reinforced, high-density polyethy¬ 
lene (HDPE) pipe that is filled with liquid scintillator. As 
shown in Fig. |T| the pipe is positioned so that it extends 
radially outward from the target station. The upstream 
location maximizes the sensitivity to oscillations by be¬ 
ing the shortest possible distance from the source, given 
spatial constraints. We have found that a long detector 
(120 m, 684 tons) is most suitable for optimizing sensitiv¬ 
ity to oscillations across a wide range of the most perti¬ 
nent parameter space, in consideration of current global 
fit results, the neutrino energy, 1/R 2 , and estimated cost. 

The interior of the pipe contains a cylinder of highly 
reflective panels, which optically separate the active vol¬ 
ume from the cosmic ray (CR) veto. Hoops of inward¬ 
facing silicon photomultipliers (SiPMs) are mounted on 
the interior of the panels. There are 100 equally-spaced 
SiPMs per hoop, and each hoop is separated longitudi¬ 
nally by 10 cm (see Fig. [2]). The space surrounding the 
inner target region on the other side of the panels is the 
10 cm-thick veto region. The surfaces of the veto region 
are painted white, or lined with a Tyvek®-like material, 
for high reflectivity. Along the innermost side of the veto 
region are 120 hoops of outward-facing SiPMs that each 
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run along the circumference of the pipe. The hoops have 
100 SiPMs each and are positioned at 1 m spacing along 
the inside of the veto region. The 10 cm spaces at the 
ends of the pipe are also instrumented. Each veto end 
cap is instrumented with 100 SiPMs that all face axially 
outward and are spaced equally apart on a circle with 
1 m radius. 

SiPMs are employed in both the target and veto re¬ 
gions because of their compact size and reduced cost 
when purchased in bulk. Currently available SiPMs typ¬ 
ically have a quantum efficiency around 30%. In order to 
further reduce cost, we plan on multiplexing the SiPM 
channels. For the target region, each channel of readout 
electronics monitors 25 out of the 100 total SiPMs on a 
hoop. For the veto region, one channel monitors one side 
or end cap hoop. The active area of a SiPM can range 
from 1 mm 2 to about 6 mm 2 . Assuming 6 mm 2 SiPMs, 
with 1200 hoops containing 100 SiPMs each, the target 
region will have a photocathode-coverage of only ~ 0.4%. 
Despite this low coverage, simulations of the experiment 
described in the next section indicate that there are an 
adequate number of SiPMs to achieve the goals of the 
experiment. 

The KPipe detector succeeds despite the sparse 
amount of instrumentation in the inner region because 
of its use of liquid scintillator as the detector medium. 
The low photocathode coverage is overcome by the large 
amount of light produced by the scintillator per unit 
of energy deposited. Scintillators under consideration 
for KPipe include those based on mineral oil and lin¬ 
ear alkylbenzene (LAB). One example of a currently- 
deployed mineral oil-based scintillator is the one used by 
the NO^A experiment m- This scintillator is a mix¬ 
ture of 95%-by-mass mineral oil with 5% pseudocumene 
(1,2,4-trimethylbenzene) along with trace amounts of 
PPO (2,5-diphenyloxazole) and bis-MSB (1,4-Bis(2- 
methylstyryl)benzene) wavelength shifters [38] . The UV 
photons emitted by the pseudocumene excite the PPO, 
which, as the primary scintillant, re-emits in the range 
of 340-380 nrn. These photons are then absorbed by the 
bis-MSB and reemitted in the 390-440 nm range. Along 
with developing their scintillator, the NOizA experiment 
has also established the methods to manufacture large 
quantities of it at a relatively low cost. Other exam¬ 
ples of mineral oil-based scintillators are those offered by 
Saint-Gobain. For reference, the light yield of these scin¬ 
tillators range from 28% to 66% of anthracene or ^4500 
to ~11400 photons/MeV [39]. Besides mineral oil, an¬ 
other option is to use a LAB-based liquid scintillator, 
similar to that being used by the SNO+ experiment m- 
This liquid scintillator consists of the LAB as solvent with 
PPO acting as the wavelength shifter. The advantage of a 
LAB-based liquid scintillator over those based on mineral 
oil is that it has a comparable light yield to the brighter 
Saint-Gobain scintillators [41] while also being less toxic. 
In order to be conservative, we assume in simulations of 
the KPipe detector (discussed in the next section) a light 
yield consistent with the dimmest mineral oil based liquid 


scintillator from Saint-Gobain (4500 photons/MeV). The 
liquid scintillator that is eventually employed for KPipe 
will be some optimization between light yield, cost, and 
safety. 


III. SIMULATION OF THE EXPERIMENTAL 
SETUP 

In order to study the performance capabilities of 
KPipe, we have created simulations of both the neutrino 
source and the detector. The source simulations, using 
both Geant4 [42] and MARS15 [43], model 3 GeV kinetic 
energy protons hitting the mercury target. The result¬ 
ing particles are propagated, and the kinematics of all the 
neutrinos produced are recorded. A semi-realistic geome¬ 
try is employed with Geant4 for the target and surround¬ 
ing material, although the majority (86%) of 236 MeV 
are found to originate within the mercury target. About 
75% of the K + are found to DAR within 25 cm of the 
upstream end of the mercury target and the ratio of 
from K + DAR to from K + decay-in-flight over 4-7T 
is ~13:1. The K + production rate varies depending on 
which simulation software is used. The Geant4 model cal¬ 
culates the 236 MeV yield to be 0.0038 v t , per proton 
on target (POT), whereas the MARS15 model predicts 
0.0072 z/^/POT. Later, when calculating the sensitivity 
of the experiment in Section [V] we will quote a sensitiv¬ 
ity which relies on the MARS15 model for kaon produc¬ 
tion, as it has been more extensively tuned to data than 
Geant4 [33] . 

The v IL flux is propagated to the KPipe detector whose 
closest end to the source is 32 m away. The z/ ;l flux for 
—0.25 < cos 6 Z < —0.16, where 6 Z is the neutrino angle 
with respect to the proton direction (+z), representative 
of the full detector length, is shown in Fig. [3] (left). The 
time distribution of all neutrinos coming from the source 
is shown in Fig. [3] (right). The two 80 ns wide proton 
pulses can be seen in the figure, while the blue histogram 
shows the neutrinos coming from kaon decay. 

The interactions of neutrinos with the detector and 
surrounding materials are modeled with the NuWro event 
generator [35], and the v„CC cross section and expected 
rate can be seen in Fig. E] Notably, the signal (KDAR) 
to background (non-KDAR) ratio is 66:1 integrated over 
all energies. In other words, if a neutrino-induced muon 
is observed, there is a 98.5% chance that it came from a 
236 MeV CC interaction. Given 5000 hours/year of 
J-PARC 1 MW operation (3.75 xlO 22 POT/year), con¬ 
sistent with Ref. [3B], we expect 1.02 x 10 5 KDAR CC 
events/year in the 684 ton active volume. 

For each generated 236 MeV CC interaction on 
carbon, NuWro provides the momentum of the outgo¬ 
ing muon and any final state nucleons (typically a single 
proton). Fig.[5]shows the kinetic energies of the resulting 
KDAR signal muons along with the non-KDAR muons. 
The v, L CC cross section on carbon at 236 MeV according 
to NuWro and employed for the event rate estimate here 
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FIG. 3: Left: The muon neutrino and antineutrino flux with —0.25 < cos 9 Z < —0.16, representative of the full detector length, 
where 9, is the neutrino angle with respect to the proton direction (+z). Right: The neutrino creation time relative to the two 
beam pulses (dotted lines). This distribution includes neutrinos emitted over all solid angles and energies. 



FIG. 4: The charged current event rate, for neutrinos with 
—0.25 < cos 9 z < —0.16, along with the employed jy, CC 
cross section. The monoenergetic 236 MeV neutrino signal 
is clearly visible above the “background” non-monoenergetic 
events, mainly coming from kaon decay-in-flight. 


is 1.3 x 10 -39 cm 2 /neutron. This is consistent with the 
Random Phase Approximation (RPA) model’s pf7| - PnT| 
cross section prediction of (1.3+0.2) x 10 -39 cm 2 /neutron 
(RPA QE+npnh). While NuWro is the only generator 
we use to produce simulated events, we did compare the 
kinematic distributions given by NuWro to that provided 
by GENIE [50] and the Martini et al. RPA model R)l . 
which includes multi-nucleon effects. We find that the 
difference in the muon kinematic predictions among the 
models is not large enough to significantly change the 
detector simulation and oscillation sensitivity results. 

Particle propagation through the detector is modeled 
using the Geant4-based simulation package RAT m- 
The detector geometry input into the simulation is as de¬ 
scribed in the previous section. The detector is assumed 
to be on the surface and is surrounded by air only. Neu¬ 
trino events are distributed over a 5 m x 5 m x 140 m box 



FIG. 5: The muon and total kinetic energy (KE to t = KE M + 
J/KEp) for the signal 236 MeV charged current events 
compared to all other ry,. Only neutrinos with —0.25 < 
cos 9 z < —0.16 are considered. The ratio of integrated sig¬ 
nal (black) to background (red) is 66:1. 


that fully contains the 120 m long, 3 m diameter cylin¬ 
drical detector. The distribution of events in the box is 
weighted to take into account the 1/R 2 dependence of 
the flux along with the density of the various materials 
in the simulation. The small divergence in the neutrino 
direction is also considered. The RAT package includes 
a model for scintillator physics that derives from models 
previously employed by other liquid scintillator experi¬ 
ments such as KamLAND. The processes that are con¬ 
sidered include scintillation, absorption, and reemission. 
All three have wavelength dependence. The reflectivity 
of surfaces in the detector is simulated using the models 
built into Geant4. 

In addition to the simulation of KDAR neutrino inter¬ 
actions with the detector and surrounding material, we 
simulate the propagation of CR throughout the volume. 
We use the simulation package CRY [52] to study the 
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CR particle flux, which generates showers consisting of 
some combination of one or more muons, pions, electrons, 
photons, neutrons, or protons. The dark rate of SiPMs 
is also included in the simulation of the SiPM response. 
We use a dark rate of 1.6 MHz for each of the 130,200 
4 mm x 4 mm SiPMs (0.4% photo-coverage) along with 
a total quantum efficiency of 30%. The dark rate comes 
from the specification for SenSL series C SiPMs which 
have an advertised dark rate of < 100, 000 Hz/mm 2 f52]. 


IV. ISOLATING AND RECONSTRUCTING u fJ 
EVENTS FROM THE KDAR SOURCE 

Signal events from the KDAR neutrino source are iden¬ 
tified by the observation of two sequential pulses of light. 
The first pulse comes from the muon and vertex en¬ 
ergy deposition. The next signal is from the Michel 
electron produced by the decay of the muon (t' /i 12 C —>■ 

X, —> e~u^i' e ). We apply a pulse finding algorithm 
to identify both light signals from the SiPMs. The algo¬ 
rithm uses a rolling 20 ns window over which the num¬ 
ber of hits in the SiPMs are summed and the expected 
dark hit contribution in the window is subtracted. The 
first pulse is found when the hit sum with subtraction 
is above a given threshold, specifically one that is four 
times larger than the standard deviation of the expected 
number of dark hits. After the first pulse is found, the 
algorithm searches for the Michel signal using the same 
method, except that the threshold is raised to account 
for both the expected dark noise and the contribution of 
SiPM hits from the first pulse. This expected hit contri¬ 
bution is dictated by the decay time of the scintillator. 
After isolating coincident signals, the position along the 
detector of both the primary interaction and Michel sig¬ 
nal is determined by the photoelectron-weighted position 
of the hits seen by the SiPMs. Using the position of the 
prompt pulse, we find that the vertex position resolution 
of the interaction is 80 cm. The current proposed readout 
is likely unable to reconstruct more detailed information 
about the event such as the muon angle, although this 
information is not necessary for KPipe’s primary mea¬ 
surement. 

Fig. § shows the number of photoelectrons (pe) in the 
first pulse as a function of total kinetic energy, KE to t, 
defined as the total kinetic energy of the muon and any 
final state protons (KE to t = KE M + E KE p ). The figure 
shows simulated data from 236 MeV KDAR CC in¬ 
teractions. The first pulse usually contains over 800 pe, 
indicating that, despite the low photocathode coverage, 
the amount of observed light for the signal events is high 
enough for efficient reconstruction. Further, the figure 
shows that KE to t correlates well with the number of pe 
seen. Using the peak of this distribution, the detector 
light yield is calculated to be 9.2 pe/ MeV, which includes 
effects from quantum efficiency, photocathode-coverage, 
and absorption. 



KE tot (MeV) 


FIG. 6: The number of photoelectrons in a 236 MeV CC 
event’s first pulse versus the total kinetic energy (KE to t = 

ke m + Eke p ). 

A. Isolating the Signal 

The primary background to the v jJt CC signal events 
comes from stopping cosmogenic muons in the detec¬ 
tor. We envision applying the following selection require¬ 
ments in order to select signal interactions and reject CR 
backgrounds: 

1. the interaction signal (prompt) occurs within 
125 ns windows following each of the two 80 ns 
beam pulses, 

2. the interaction signal has a reconstructed energy in 
the range 22 < E v i s < 142 MeV (200 < pe < 1300), 

3. the Michel signal occurs within 10 p,s of the first 
pulse, 

4. the Michel signal reconstructed visible energy is 
11 < E vis < 82 MeV (100 < pe < 750), 

5. the distance between the interaction signal and the 
Michel signal is less than 1.5 m, and 

6. the summed pulse height in the ten nearest veto 
SiPM hoops to the interaction signal is less than 
four times the dark rate a within a 125 ns window 
after the start of each 80 ns beam pulse. 

Note that for the cuts on visible energy, E v j s , the corre¬ 
sponding values in pe are given in parentheses. These are 
the values used in the Monte Carlo study of the KDAR 
signal efficiency and CR background rejection. 

The first cut (1) takes advantage of the pulsed proton 
beam. Accepting events only within a 125 ns window 
after each 80 ns proton pulse efficiently selects 99.9% of 
the KDAR neutrinos while removing many of the events 
coming from other neutrino sources. The small 125 ns 
event window also limits the rate of CR ray events even 
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FIG. 7: The total kinetic energy (KE to t = KE,, + ^KE P ) 
versus the energy of neutrinos from CC interactions in KPipe. 
Only neutrinos with —0.25 < cos# z < —0.16 are considered. 
The Z-axis units are arbitrary. 


before the other selection cuts are applied. According to 
the simulation, CR particles create at least one detectable 
flash in either the target region or veto in only 0.87% of 
all windows. 

The second cut (2) utilizes the fact that, because the 
signal events come from monoenergetic neutrinos, the 
energy of the outgoing particles falls in a fairly narrow 
range. Fig. [7] shows the total kinetic energy of the muons 
and any final state protons, KE to tj as a function of neu¬ 
trino energy for CC events in the detector. The up¬ 
per bound of 142 MeV ensures that the signal neutrino 
events are preserved with high efficiency, while remov¬ 
ing non-KDAR muon neutrinos at higher energies. More 
importantly, the upper bound removes bright CR events. 
Based on the simulation, 72% of all detectable CR events 
(i.e. ones that produce one or more detected flashes) 
are removed by the high energy cut, many of which are 
through-going muons. Along with kaon decay-in-flight 
neutrinos, the low energy bound also removes all rele¬ 
vant backgrounds from CR-induced spallation products 
and is well above the visible energy from radiogenic back¬ 
grounds. With both a high and low energy cut on the 
first pulse, 87% of all CR events are removed. 

The cuts related to Michel electron timing, energy, and 
spatial coincidence (cuts 3-5) are chosen to efficiently re¬ 
tain signal while removing most of the in-time through- 
going CR muons that traverse the detector, as well as 
other backgrounds. A coincident signal coming from non¬ 
stopping muons can occur due to a CR shower with two or 
more particles or an associated muon spallation-induced 
isotope. The timing, energy, and spatial cuts on the 
Michel candidate reduce much of this coincident back¬ 
ground. Applying the above cuts along with the Michel 
pulse cuts reduces the CR rate to 750 Hz, which means 
that only 0.01% of all signal windows will contain a CR 
event. At this stage in the cuts, less than two percent of 
detectable CR events remain. 


The final cut (6) applied removes all events that create 
a flash of light in the veto. The veto is only 10 cm thick 
and is more sparsely instrumented than the target region. 
However, enough light is produced that the veto is able 
to reject 99.5% of all detectable CR events with at least 
one muon. We find that lining the walls of the veto with a 
highly reflective material plays an important role in the 
veto performance. With all cuts applied, we estimate 
that the rate of CR events is 27 Hz over the entire active 
volume. 

In addition to CR backgrounds and non-KDAR muon 
neutrino events, an additional coincident background can 
come from beam-induced neutron interactions that pro¬ 
duce a A + in the detector that subsequently decays into 
a 7r + . The latter can then stop and decay to a muon fol¬ 
lowed by a Michel electron. We assume that this back¬ 
ground is negligible for this study. All in-time beam- 
related backgrounds will be measured before deploying 
KPipe, and adequate shielding will be installed in order 
to mitigate them. 

Overall, our studies indicate that the dominant back¬ 
ground is from CR shower events that are not removed 
by the above cuts. Of the 27 Hz rate that passes, the 
simulations show that 70% of the rate is due to stopping 
muons. The remaining 30% is due to showers involving 
photons, electrons, and neutrons. In the simulation, we 
do not include any additional passive shielding, for ex¬ 
ample coming from overburden. If the detector is buried 
or shielded, we expect these non-muon backgrounds to 
be further reduced. The CR background should be dis¬ 
tributed uniformly throughout the detector and can be 
measured precisely using identified out-of-time stopped 
muons. As a result, only the statistical error from the 
total number of background events expected to pass the 
cuts is included in the sensitivity analysis, described later 
in Section El 


B. Detection efficiency 

The cuts introduce inefficiency in the signal. We as¬ 
sume that the neutrino events are distributed evenly in 
radius and fall as 1 /R 2 throughout the detector. Signal 
events near the lateral edge of the target region can exit 
the detector before the muon can decay. This leads to 
an acceptance that is a function of radius. Based on an 
active detector radius of 1.45 m, we find an acceptance 
of 87% with respect to KDAR v^ CC interactions whose 
true vertex is in the target region. The selection cuts 
described above are 89% efficient according to the simu¬ 
lation. This includes events where the muon is captured 
by the nucleus, which occurs in the target region 6% of 
the time. For a subset of these events, there is also an 
additional 0.75% dead-time loss due to the rate of CR 
events in the veto. 

In summary, the total efficiency for all signal events 
is 77%, leading to an expected total KDAR CC rate 
of 7.8 x 10 4 events distributed along the pipe’s active 
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volume per year of running. This is on average 4.9 x 
10 ~ 5 KDAR events per proton beam window without 
oscillations. This compares with 3.4 x 1CU 6 CR events per 
proton beam window. In the most upstream 1 m of the 
detector, the unoscillated signal to background ratio is 
about 60:1; in the most downstream 1 m of the detector, 
the unoscillated signal to background ratio is about 3:1. 


V. SENSITIVITY 

The oscillation probabilities for three different Am 2 
values (1, 5, 20 eV 2 ) can be seen in Fig [ 8 j The error bars 
correspond to the statistical uncertainty associated with 
a 3 year measurement with a CR rate of 27 Hz. This 
background rate corresponds to 132 CR events that pass 
our selection cuts for each 1 m slice of the detector. 

The sensitivity of the experiment is evaluated using 
a shape-only y 2 statistic similar to that described in 
Ref. [54]. However, we replace the covariance matrix with 
the Neyman y 2 convention, since we do not include any 
correlated systematic uncertainties between each L/E 
bin. Using Eq. [ 5 ] for the oscillation probability, the y 2 
value at each pair of oscillation parameters, Am 2 and 

4 , is calculated by comparing the no-oscillation signal 
(N/’ un + N dkgd ) to the oscillation signal (Np osc + N dkgd ) 
in each L/E bin, i. Here, N/' un and ND 0 ^ are defined 
as the number of expected events in bin i given a no¬ 
oscillation prediction and an oscillation prediction, re¬ 
spectively. The number of events in a bin due to back¬ 
ground is then added to the prediction. The A L value 
used in setting the bin size is 0.5 m. Defining for each 
i th L/E bin the difference between the no-oscillation and 
oscillation signal, n,, where 

n; = (ND un + N dkgd ) - (£ND osc + N dkgd ) , (4) 

the y 2 is then 

nbins 2 

x 2 = y — 111 .. . . (5) 

^ j^-bkgd v ' 

The normalization constant, £, in Eq. [4] is included in 
order to make the analysis shape-only and is constrained 
to be 


£ 


e N r' 


E N r 


( 6 ) 


For the 90% confidence limit reported, a one degree of 
freedom, one-sided raster scan threshold of y 2 =1.64 is 
used. The 5cr threshold is y 2 =25.0, considering a one 
degree of freedom, two-sided raster scan. 

For the subsequent sensitivity plots, the oscillation pre¬ 
diction, N/’ osc , has been simplified by the two flavor ap¬ 
proximation to the 3+1 neutrino oscillation model (Equa¬ 
tion]^]), where we define sin 2 ( 20 M/i ) = 4 1 C/^, 41 2 (1 — 11 /^ 41 2 )- 



FIG. 8: Three sample oscillation probability measurements 
as a function of L/E for 3 years of running. The error bars 
incorporate statistical uncertainties of both the iy, signal and 
the cosmic ray background. 


Parameter 

Value 

Detector length 

120 m 

Active detector radius 

1.45 m 

Closest distance to source 

32 m 

Liquid scintillator density 

0.863 g/cm 3 

Active detector mass 

684 tons 

Proton rate (1 MW) 

3.75 xlO 22 POT/year 

KDAR yield (MARS15) 

0.0072 i^/POT 

iz,, CC a @ 236 MeV (NuWro) 

1.3 x 10~ 39 cm 2 /neutron 

Raw KDAR CC event rate 

1.02 x 10 5 events/year 

KDAR signal efficiency 

77% 

Vertex resolution 

80 cm 

Light yield 

4500 photons/MeV 

i/^ creation point uncertainty 

25 cm 

Cosmic ray background rate 

27 Hz 


TABLE I: Summary of the relevant experimental parameters. 


The KPipe search for sterile neutrinos, which uses only 
the relative rate of events along the pipe, is helped by the 
fact that uncertainties associated with the absolute nor¬ 
malization of the event rate expectation are not relevant 
for this shape-only analysis. This includes theoretical 
uncertainties in the kaon production and neutrino cross 
section. Instead, the dominant uncertainty associated 
with the weight of each bin comes from the combined 
statistical uncertainty of the u fl measurement and the 
CR background. In the sensitivity studies, we assume a 
CR background rate of 27 Hz over the entire detector. 
Further, there are two uncertainties associated with the 
neutrino baseline L: the creation point of the v )t from 
the decaying K + has an uncertainty of 25 cm; the recon¬ 
structed position resolution, described in Section [IV] has 
an uncertainty of 80 cm. There is no uncertainty asso- 




















































FIG. 9: The projected sensitivity of KPipe to muon neutrino 
disappearance with 3 years of running, including the cosmic 
ray background, signal efficiencies, and reconstruction uncer¬ 
tainties described in the text. The red contours are the global 
allowed regions given by Collin et al. [55]. 



FIG. 10: The 90% CL sensitivity of KPipe with 6 years of 
running, compared to the sensitivity from 6 years of the SBN 
program. The KPipe sensitivity estimate includes the cosmic 
ray background, signal efficiencies, and reconstruction uncer¬ 
tainties described in the text. 


ciated with the energy reconstruction since the u,, have 
a definite energy. We also include a total detection effi¬ 
ciency due to the selection cuts, d ead-time, and escaping 
muons described in Section IV A of 77%. A summary of 
the relevant experimental parameters and assumptions 
can be seen in Table Q] 


Fig- m shows the projected 90% and 5(7 sensitivity of 
KPipe to v fJ: —>• for 3 years of running. The global fit 
allowed regions, given in red, were produced using a new 
software package based on the previous work of Ignarra 
et al. [IT] . We refer to this work as “Collin et al .” [55] , 
The fit includes the datasets described in Ref. [5B] with 
the exception of the atmospheric limit. The model pa¬ 
rameters are explored using a Markov chain Monte-Carlo 
algorithm. Contours are drawn in a two-dimensional pa¬ 
rameter space using 2 degree of freedom % 2 values for 
90% and 99% probability. After 3 years of KPipe run¬ 
ning, the 5 ct exclusion contour covers the best fit point 
at Am 2 = 0.93 eV 2 and sin 2 (26 ) jLtAl ) = 0.11. 

Fig. [TO] shows a comparison between KPipe’s predicted 
six year 90% sensitivity and the predicted sensitivity of 
SBN |22l assuming 6.6 x 10 20 POT (3 years) in SBND 
and the ICARUS-T600 and 13.2 x 10 20 POT (6 years) in 
MicroBooNE. The dashed contour represents the com¬ 
bined 90% excluded region based on the muon neutrino 
disappearance results of MiniBooNE and SciBooNE m 
SBN and KPipe have similar sensitivity reach in the 
Am 2 = 1 — 4 eV 2 region, however SBN performs better at 
low-Am 2 and KPipe at high-Am 2 ; the complementarity 
between the experiments is clear. 


VI. CONCLUSION 

The J-PARC MLF facility provides a unique and in¬ 
tense source of neutrinos in the form of monoenergetic 
236 MeV muon neutrinos coming from the decay-at-rest 
of positively charged kaons. The KPipe experiment seeks 
to take advantage of this source for a decisive v tJ dis¬ 
appearance search at high-Am 2 in order to address the 
existing anomalies in this parameter space. The 120 m 
long, 3 m diameter liquid scintillator based active volume 
(684 ton) will feature 0.4% photo-coverage for detecting 
these Vfi CC events in an attempt to discern an oscillation 
wave along the length of the detector. 

In contrast to other neutrino sources, the KPipe neutri¬ 
nos are dominantly monoenergetic. This provides a great 
advantage in searching for neutrino oscillations. A neu¬ 
trino (or antineutrino) induced double-coincidence muon 
signal detected with KPipe has a 98.5% chance of being 
from a 236 MeV CC event. This simple fact allows 
the active detector requirements to be extremely modest, 
the systematic uncertainties to be practically eliminated, 
and the detector’s energy resolution to be only a weak 
consideration. 

Within three years of running, KPipe will be able to 
cover the current global fit allowed region to 5(7. The 
sensitivity for a 6 year run at the J-PARC facility will 
enhance existing single experiment limits on disap¬ 
pearance by an order of magnitude in Am 2 . Such a mea¬ 
surement, when considered alone, or in combination with 
existing and proposed electron flavor disappearance and 
appearance measurements, can severely constrain mod¬ 
els associated with oscillations involving one or more light 
sterile neutrinos. 
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